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Using an integration formula recently derived by Conrey, Farmer and Zirnbauer, we calculate the 
expectation value of the phase factor of the fermion determinant for the staggered lattice QCD action 
in one dimension. We show that the chemical potential can be absorbed into the quark masses; the 
theory is in the same chiral symmetry class as QCD in three dimensions at zero chemical potential. 
' In the limit of a large number of colors and fixed number of lattice points, chiral symmetry is broken 

spontaneously, and our results are in agreement with expressions based on a chiral Lagrangian. In 
this limit, the eigenvalues of the Dirac operator are correlated according to random matrix theory 
£S) ' for QCD in three dimensions. The discontinuity of the chiral condensate is due to an alternative 

to the Banks-Casher formula recently discovered for QCD in four dimensions at nonzero chemical 
potential. The effect of temperature on the average phase factor is discussed in a schematic random 
matrix model. 
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I. INTRODUCTION 



In spite of significant recent progress, QCD at nonzero chemical potential remains a notoriously hard problem 
[USEl- In particular, first principle nonperturbative results at low temperature are absent because the phase of the 
fermion determinant invalidates probabilistic methods to evaluate the partition function. This problem is known as 
the sign problem. The sign problem is particularly severe if the phase of the fermion determinant results in a different 
free energy, i.e. 

(N - 

(Idet^DI) f Pq 

If the phase quenched free energy, F pq differs from the free energy of QCD with Nf flavors, Fn,, the number of 
required gauge field configurations grows exponentially with the volume. The ratio defined in equation (|TJ) can also 
, be interpreted as the average phase factor of the fermion determinant 0, [1[ 

O : {p iN f es (e m f e \det N W\) 

{ )pq ~ (\det N m\) ■ () 

Alternatively, one can define the average phase factor with respect to the full QCD partition function [fj, [f| 

(e )Nt= (det^D) ' (3) 

which has also been used in the literature [1, [E fl & ll> S El EL G2] as a measure for the severity of the sign problem. 

When the chemical potential is sufficiently small, the average phase factor can be studied by means of chiral 
perturbation theory 0,0. Recently, the average phase factor was analyzed in the microscopic domain of QCD (l3l.[l3|. 
where only the constant fields in the chiral Lagrangian contribute to the mass and chemical potential dependence of 
the partition function. It was found that the sign problem is not serious for /i < m^/l- For fi > m x /2 the chiral 
condensate of the phase quenched theory rotates into a pion condensate resulting in a free energy that is different 
from the full theory and a severe sign problem. 

In this paper we study the sign problem for Euclidean QCD in one dimension. In one dimension, the only effect of 
the gauge field is in the boundary conditions, and lattice QCD simplifies to a matrix integral. A general formula for 
exactly this type of integrals was recently derived by Conrey, Farmer and Zirnbauer [lj| . Using what we will call the 
CFZ formula, exact analytical results for the one dimensional QCD partition function and the average phase factor 
will be obtained. 

The study of QCD in one dimension at finite chemical potential has had a long histor y. O ne reason to study this 
model is that it is an effective model for the large N c limit of strong coupling QCD [Ttj Il7l. [Til . [l9l . I20I |2lj . Among 
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others, it has been successfully used to explain [z| HU HH puzzling lattice results for quenched QCD at nonzero 
chemical potential. More recently, one dimensional QCD was used to study the complex zeros of the partition 
function 0,[2i, [2^ . 

The sign problem of QCD at nonzero chemical potential arises because of the nonhcrmiticity of the Dirac operator. 
QCD in one dimension is not truly nonhermitean: Instead of being scattered in a two-dimensional domain of the 
complex plane [13], the eigenvalues of the Dirac operator are localized on an ellipse in the complex plane [H, |22T |. 
Another manifestation of the mild nonhermiticity is that the chemical potential can be absorbed into the quark 
masses. For each flavor in one dimension with quark mass m, we can associate flavors with mass m + fj, and — m + fi. 
In the large N c limit, this results in spontaneous chiral symmetry breaking according to U(2Nt) — ► U(Nf) x U(Nf). 
Therefore, QCD in one dimension is in the same chiral symmetry class as QCD in three dimensions. The microscopic 
limit of QCD in one dimension is equivalent to the microscopic limit of QCD in three dimensions [28j |. Another 
peculiarity of QCD in one dimension is that, already in the free theory, the eigenvalues near zero are spaced inversely 
proportional to the volume. Therefore, in the thermodynamic limit (i.e. at zero temperature) the chiral condensate 
is discontinuous across the imaginary axis. However, this type of symmetry breaking is not a collective phenomenon, 
and there are no associated Goldstone bosons. It is reminiscent to the alternative to Goldstone's theorem proposed 
by McKane and Stone [29) }. Spontaneous symmetry breaking with Goldstone bosons takes place for a fixed number 
of lattice points in the limit N c — > 00. 

A chemical potential excites color singlet excitations with baryon charge qu and mass when q^fi > Mk- This is 
also the case for one-dimensional QCD. For gauge group U(N C ), the only color singlet excitations are mesons that are 
uncharged with respect to the chemical potential, whereas for SU (N c ) gauge group we have both neutral mesonic and 
charged baryonic color sin glet excitations. The complex conjugate of the fermion determinant can be interpreted in 
terms of conjugate quarks |23l. [30j which have a baryon charge that is opposite to that of regular quarks. Therefore, 
both for U(N C ) and SU{N C ), the phase quenched partition function, where the fermion determinant has been replaced 
by its absolute value, has charged mesonic excitations made out of quarks and conjugate anti-quarks. This results in 
a phase transition at /z = \i c = The SU(N C ) theory with Nf flavors has only charged baryonic excitations and 

will have a phase transition at [i = fj, c = m,B/N c . In one dimension it turns out that the pion and the baryon have 
the same mass per quark number so that the critical chemical potential of the phase quenched U(N C ) theory and the 
SU(N C ) theory is the same. 

For the U(N C ) partition function we expect the sign problem to be severe when fi > fj, c because the phase quenched 
partition function has a phase transition at fi = /i c whereas the normal theory remains in the same phase. For 
SU(N C ), both the normal partition function and the phase quenched partition function do have a phase transition to 
a phase of free quarks at fj, = [i c and have the same free energy not only for fi < [i c but also for /i > /i c . Therefore, 
in one dimension there is no severe sign problem for SU(N C ) QCD. 

Because one dimensional QCD with gauge group U(N C ) does not have a phase transition, the chiral condensate 
in the thermodynamic limit or large N c limit is discontinuous at m = 0, independent of the value of the chemical 
potential. Because the Dirac eigenvalues are located on an ellipse for /j, ^ 0, this discontinuity cannot be related to 
the Dirac spectrum by means of the Banks-Casher formula. In [31] a different mechanism to explain the discontinuity 
was discovered . It was found that an oscillating contribution to the spectral density with an amplitude that diverges 
exponentially with the volume, is responsible for the discontinuity of the chiral condensate. We will show that a 
similar mechanism is at work for U (N c ) gauge theory in one dimension. 

Lattice QCD in one dimension will be introduced in section II, where we also discuss its continuum limit, mean field 
results and the Conrey-Farmer-Zirnbauer formula. In section III we will evaluate the U(N C ) partition function and 
its phase quenched version. The average phase factor for U (N c ) is calculated in section IV. In section V we study one 
dimensional QCD with SU(N C ) as gauge group and evaluate the regular and phase quenched partition functions and 
the average phase factor. The effect of temperature will be illustrated with results from a schematic random matrix 
model in section VI. The connection with the microscopic domain of QCD in three dimensions is discussed in section 
VII, and the Dirac spectrum and its relation with chiral symmetry breaking is analyzed in section VIII. Concluding 
remarks are made in section IX. 
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II. QCD IN ONE DIMENSION 
A. Lattice QCD 



The staggered lattice QCD Dirac operator in one dimension is given by 



D = 



ml 



-e-v-Uyi ml 



-"14/2 \ 




ml 

e-»U*_ ln /2 



e»U n . ln /2 
ml 



(4) 



/ 



We have used anti-periodic boundary conditions and the gauge helds on the links are taken to be in U(N C ) or SU (N c 
The chemical potential /j, is an imaginary vector field introduced according to the Hasenfratz-Karsch prescription [32 
In general, the lattice QCD partition function for Nt flavors with mass m is given by 



(5) 



where the product is over the links and Sym is the Yang-Mills action. In one dimension, the partition function 
simplifies substantially. First, there is no Yang-Mills action, and second, by unitary transformations the integral over 
the links can be reduced to a single integral. In the gauge where all gauge fields except U n i = U are equal the unity, 
the fermion determinant reduces to [24| 



dct D = 2 



— nN c 



det[e 



-nuc _|_ e n ^JJ . 



with fi c given by 



fx c = sinh 1 ■ 



(6) 



(7) 



This value will turn out to be the critical value of the chemical potential. From now on, the overall factor 2~ nNa will 
be absorbed into the normalization of the partition function. To avoid sign factors, n is taken to be even throughout 
this paper. 

The standard method to evaluate the QCD partition function in one dimension is to use the eigenvalues of U as 
integration variables 0, |b(l [13, d, HH . For example, one finds this way that the result for Nf — 1 with U(N C ) as 
gauge group is given by [24j | 



Zjv /= i(/i c , jit) = / dUdctD = 

Ju{N c ) 



sinh((^ e + l)nn c ) 
sinh(n/Lt c ) 



(8) 



For most partition functions that are considered in this paper it is not possible to obtain analytical results by means 
of this method. Instead we will use powerful integration formulas for unitary integrals that were recently derived in 
(TH - These integrals are based on the color-flavor transformation [33[ which has also been applied to lattice QCD 
with baryons in the canonical ensemble [H, [35| . 



B. Continuum Theory 



The continuum limit of the staggered lattice action is given by 

n _ ( m d + iA + fi \ , , 

^co„t ~ y do + lAo + fI m J ' W 

where Ao is a Hermitean N c x N c matrix and the off-diagonal blocks connect even and odd lattice sites. What is 
special in one dimension is that the off-diagonal blocks are identical. The eigenvalue equation (do + iAojipk = iuk^pk 
is solved by 



$ = Pe- { J? A « dt Xo with d oXQ = iE X o, 



(10) 
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with Pe the path ordered exponent. Nontrivial eigenvalues are obtained by imposing boundary conditions on ip. At 
nonzero chemical potential the eigenvalues of the Dirac operator are given by 

Afc = to ± (iuk + p) with Uk £ M. (11) 

This is very different from QCD with d > 2 where the eigenvalues of the QCD Dirac operator at ^ ^ are scattered in 
the complex plane. Also the eigenvalues of the one dimensional staggered lattice Dirac at nonzero chemical potential 
are localized on a curve in the complex plane. 

Another consequence of the structure of the Dirac operator (|9|) is that the fermion determinant can be rewritten as 

det Dcont = det(<9o + + M + m ) det(<9o + iA + fi — m), (12) 

which can be interpreted as a two-flavor partition function with masses m + fi and — m + /i. This Dirac operator 
has the same structure as QCD in three dimensions: In the large N c limit, chiral symmetry is broken spontaneously 
according to £7(2) — > U(l) x U(l) with the squared Goldstone masses proportional to the difference of the positive 
and negative quark masses, i.e. ~ m — (—to) = 2m. At low energy, these Goldstone modes interact according to a 
chiral Lagrangian determined by the pattern of chiral symmetry breaking. In the large N c limit and m^P -C 1 (with 
(3 the length of the one-dimensional box) the QCD partition function in one dimension is therefore equivalent to the 
low-energy limit of QCD in three dimensions. In section [VIII this will be worked out explicitly for the microscopic 
limit of the partition function. This is the limit 

p(0) — » oo with mnp(O) = fixed and pirp(0) = fixed, (13) 

where p(0) is density of the eigenvalues, or their projections onto the imaginary axis, close to zero. Since p(0) = nN c /2w 
(see section VIII), this is the limit nN c — > oo with nN c m and nN c p fixed. When we use the term microscopic limit, 
we always mean the universal microscopic limit. This is the limit (|13[) associated with the formation of Goldstone 
bosons, i.e. the limit N c — > oo at fixed n. 



C. Mean Field Limit for Large N c 



In the large N c limit, chiral symmetry is broken spontaneously even in one dimension, so that its low-energy limit 
is a theory of Goldstone bosons. In this section we give general arguments that determine the chemical potential 
dependence of the partition function in the microscopic domain where the Compton wavelength of the Goldstone 
bosons is much larger than j3. 

As was argued in 0, Q , in the microscopic domain, the mean field limit of the partition function is given by 

Z =;(n^)«-, ("> 

where J is the value of the integration measure at the saddle point, F is the free energy, V is the space time volume, 
and ttItt are the masses of the Goldstone bosons. Let us apply this result to the average phase factor in one dimension. 
For n < n c , in the limit T — > 0, only the vacuum state contributes to the partition function, so that the free energy 
is independent of /i. For quark mass m, the equivalent QCD3 mass matrix of the Dirac operator D in (|12p is given 
by diag(— m + fi, m + /x), whereas the hermitean conjugate Dirac operator D\ has diag(— m — fjL,m — p) as equivalent 
QCD3 mass matrix. For Nf flavors, the average phase factor is the partition function with Nf + 1 fermionic quarks 
and one conjugate bosonic quark. We thus have 2(Nf + l) 2 Goldstone bosons made out of two fermionic quarks with 
squared mass m\ = 2mG (with G a constant), &{Nf + 1) fermionic Goldstone modes, half of them with squared mass 
m\ = 2(m — p)G, and the other half with squared mass m\ — 2(m + p)G. Finally, we have 2 Goldstone bosons 
composed out of two bosonic quarks with squared mass equal to m\ — 2mG. The fermionic partition function Zff f 
has 2Nj Goldstone bosons all with squared mass to 2 . = 2mG. Using (|14[) we thus find 

{e z Nfi p c ,») = { 1 -^) for ^ (15) 

For phase quenched partition functions with Nf flavors and Nf conjugate flavors, the equivalent QCD3 symmetry 
breaking pattern is U(4Nf)/(U(2Nf) x U(2Nf)), so that the total number of Goldstone bosons is equal to 8./V 2 . Of 
these, 4iV? Goldstone bosons have squared mass to 2 = 2toG, half made out of two quarks and the other half out of two 
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conjugate quarks. The other ANj Goldstone bosons are composed out of a quark and a conjugate anti-quark quark, 
half with squared mass = 2 (to — /i)G and the other half with squared mass = 2(m + /x)G. All equivalent QCD3 
Goldstone bosons of the full QCD partition function with 2Nf flavors have a squared mass equal to = 2mG. For 
fi < /i c the free energy of the phase quenched partition function and the full QCD partition function is the same so 
that, using (fT4|) . the phase quenched average phase factor is given by 



/ 2i6\ Z N +1+(N1 y([l c ,[l) t n 2 \ / 

Below we will show that the results derived in this section also follow from the zero temperature microscopic limit 
of the exact evaluation of the partition function. 

D. The Conrey-Farmer-Zirnbauer formula 

Exactly the integrals that are required for the evaluation of the average phase factor of QCD at nonzero chemical 
potential in one dimension were studied in a recent paper by Conrey, Farmer and Zirnbauer [151 ] . They considered 
the partition function 

vtr > f Jrr Adet(l-e*C/) ^j 9 det(l-e-^[/t) 

Ju(N c ) fj{ det(l - e*'U) det(l - e 4>'W) 

with dU the Haar measure of U(N C ) and ip k , <t>k complex parameters with Re(<f>j) < < Re(^). Using the color 
flavor transformation [33 | and Howe's theory of supersymmetric dual pairs, they derived the following formula 



p+q 



K£S p+q /(S p xS q ) l+l 3- 



The sum is over permutations in S p + q /(S p x S q ) that interchange any of the ip±, • • • , ip p with any of the ip p +i, ■ ■ ■ , 4> P + q - 
Partition functions with an unequal number of bosonic and fermionic determinants can be obtained from special 

limits of (|18|) . In the case of only fermionic determinants, an equivalent expression was first obtained in (3q . [37| . 

Expressions for degenerate parameters can be derived by carefully taking limits of the above formula. 

Unitary matrix integrals can also be calculated by using an eigenvalue representation of the unitary matrices. The 

integrals we are interested in are of the form 

Z= [ dUY[F(e ie »), (19) 

JU(N a ) k 

where exp(i9k) are the eigenvalues of U. Using an eigenvalue representation of the unitary matrices, the orthogonal 
polynomial method can be used to express Z as [l7| 

Z = det(B k _ l ) k<l=0 ,-,N c -i, (20) 

where 

B k = ±- r d6e tk9 F({e t9 }). (21) 
2tt J -tv 

This is the method that was used in the literature on one-dimensional QCD prior to this paper 0, [l(| [13, [H, HH • 
In a few cases the determinant in (|20|) could be evaluated explicitly resulting in expressions that are similar to those 
derived directly from the CFZ-formula. We have used (f2"0"l) to numerically check the results obtained by means of the 
CFZ-formula. 

III. EXACT EVALUATION OF THE ONE-DIMENSIONAL U(N C ) PARTITION FUNCTION 

A. Partition Function for arbitrary Nf 



In this section we evaluate the one-dimensional U(N C ) QCD partition function for Nf flavors. 



G 



To apply the CFZ formula (fTSf we rewrite the determinant (J6j> as 

det£> = e^ nN «&et{l + e n ^-^U) det(l + e^ - " - "^*/*). 
For arbitrary Nf we then find the remarkably simple answer 

e -/V c m„ (+ fc) 



Z Nf (ii c ,n)= / <fl7det JV >I> = 

JU(N C ) 



E 



nn r 

/c=i z=i 



exp(m (T (_ ;) -m CT ( +fc) ) 



with 



m_ fe 



-nfick 



and 



(22) 



(23) 



(24) 



The sum is over all permutations that interchange positive and negative masses. Notice that the /^-dependence has 
canceled from this expression. This also follows from an expansion of the determinant in powers of U and U' . Only 
terms with an equal number of factors U and factors U> are non- vanishing. 

Although the result for degenerate positive and negative masses can be obtained by carefully taking limits of 
it is simpler to start from a different representation [3a , |37| of partition function (123|) given by 



III 



<k<l<2N f 



(e M < 



1 e 1 "- 1 



1 e m - N f 
1 e m+1 



,{N f -l)m-i JN c +Nf)m-i 



,(N a +2N f -l)<m,- ! 



e (A r / -l)m_ JV/ (N e +Nj)m-N f ... e (JV c +2A r J -l)m_ JV/ 
e (JV / -l)m +1 e (JV c +JV / )m + i ... e {N e +2N f -X)m+i 



1 e ra +»/ ... e i N f~ l ) m +N f e (N c +N f )m +Nf _ e (Ar c+ 2AT / -l) m+N/ 



(25) 



For convenience we have introduced the mass matrix = (mi i, • • • , m- jv, , ni+i, • • • , m + jv, ). To obtain an expres- 
sion for degenerate masses, we Taylor expand the exponential functions exp(m + fc) and exp(m_fc) to order Nf — 1 
about exp(— nfi c ) and exp(n/x c ), respectively, and write the resulting matrix as the product of two matrices, one 
containing the Taylor coefficients, and the other containing powers of the expansion parameters. The determinant of 
the second matrix can be written as a Vandermonde determinant which cancels against part of the prefactor in (I25p . 
Our final expression for the partition function with degenerate masses is given by 

i 1 



e —2nfi c 

e -n{Nf-X)ii c 
e -n{N c +N f )fi c 



//-l e -n(JV r l) Mc 
g N J- 1 e -n(N c +N f )^ c 



g2n/i c 

e n(N f -l)i* c 
e n(N c +N f )n c 



§ Nf-l e n(Nf-l),M c 
( J J V- 1 e n(JV +AT / )/ ic 



-n(N c +2N f -l)Vc ... fi N f- l e -n(N c +2N s -l) l t a e n(N a +2N s -l)tt e ... e n(N c +2N f -l)|» c 



(26) 



with 



<5_ 



d(-nn c ) ' 



(27) 



The fe'th column (with k < Nf) is given by the S 1 ^ 1 ' derivative of the first column, and the Nf + k'th column is 
given by the 5+ ^ derivative of the TV^'th column. 

From this result one can easily derive explicit expressions for small values of Nf. The partition function for Nf = 2 

( e n(JV c +2) Mc _ e -n(N c +2)fi c \2 (N c + 2) 2 



reads 



e — nij, c \4, 



(e 1 



(28) 
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The microscopic limit of this partition function is given by 



^° 2 (McM) = [ ~ ^3 -■ (29) 

For Nf = 3 the partition function is given by: 

I n(N c +3)iJ, c _ p -n(N c +3)/j, c \3 ( n(N c + l)^ c _ p n(N c + l)^ c \ 

^/=»0*«M) = { (e n Me - e -nM,)9 J ~ + ^ (e » Me _ e -nM e )7 " ^ 

(- e n(Af c +3)Mc _ g-n(.iVo+3)/*c\ J f n(JV c +3)p c _ -n(JV c +3)Mc) 
+ 2 ) " (e nM.- e -~M.)7 " - 4 + 2 ) 2 ^ + V - (e ^_ e -^)5 

(31) 

B. The Phase Quenched Partition Function 

To calculate the average phase factor according to the definition Q for Nf — 2, we also need the one-dimensional 
phase quenched QCD partition function for two flavors which will be evaluated in this subsection. It is defined by 

Z 1+1 .(/i c ,/i) = / dU del D det D ] (32) 

Ju&U(N c ) 

= e^ inN " [ d[/det(l-e™^-^)[/)det(l-e™ ( -' I -^ ) C/ t )det(l-e™ ( ' I -^ ) ?7 t )det(l-e" ( -^^ ) f/). 

This partition function is easily evaluated using the integration formulae of [l5j]. We find 

cosh(2n(iV c + 2) He) cosh(2n(7V c + 2)/i)) 



8 sinh(n(^ c — /i)) sinh(n(/i + /i c )) sinh (n/i c ) 8 smh(n(fi — fi c )) sinh(n(/H + fi c )) sinh (n/i) 
f 



8 sinh 2 (n/i) sinh 2 (n/x c ) 

For /x = this result agrees with the Nf = 2 partition function given in Q28[) . Its microscopic limit simplifies to 



(33) 



ZT+i°(^c, M) = 16n 4 (Ai 2 _ m 2 )ai 2 + i6nVMl ~ 16n 4 (/i2 - ^2 ■ ( 34 ) 

In section IVIII we will show that this result is equal to the microscopic limit of the QCD3 partition function with 
masses —(p c + fi), -((J, c - m)j Mc - Mi Mc + 

The large iV c limit of the phase quenched partition function is given by 

e 2n(JV c +2)p c 

Zi + i*(fJ, c ,fl) = ^33-7777 -^3Tr7T77 I .. vT73rr7T7 \\2 for M < Mc, 



16 sinh(n(/x c — /i)) sinh(n(/x + /i c ))(sinh(n/i c )) 5 

e 2n(7V c +2) M 

16 sinh(n(^ — /z c )) smh(n(/x + /i c ))(sinh(n/i)) 2 

TV e 2n(Af c +2)Mc 



Z 1+1 *(n c ,n) = 7S3377--r- — .. ^ztttt^ 1 .. xTT^rrzrT^a for m > Mc 



Zi +1 ,(Mc,/x) = — -—77 -[1 + OiN- 1 )} for M = /x c , (35) 

8 sum (n/i c ) sinn(zn/x c J 

which will be used to calculate the phase quenched average phase factor in this limit. 



IV. AVERAGE PHASE FACTOR FOR U(N C ) 



In this section we will evaluate the average phase factor, first from the ratio of the full QCD partition function and 
the phase quenched partition function and in the next subsection starting from the definition in Eq. ([3]). 



8 



A. Phase Quenched Average Phase Factor 

For Nf — 2, the average phase factor with the absolute value of the fermion determinant as weight is given by 



2i8 



e- 



(e 2 ' 9 |detQD)| 2 ) Z Nf=2 (» c ,fi) 
(|det(I>)P) Z 1+1 .0*c,M)* 



(36) 



The two-flavor partition function was evaluated in subsection III A, and the phase quenched partition function was 
calculated in subsection III B. The large iV c limit of the phase quenched average phase factor is given by 

sinh i (nAIc) tor fx < [i c , 

(e 2ie ) VQ = I ^h 2 (» M ) sinh(n(^- Mc )) S inh(n( M+Pc )) 2n(jV c +2) ( M ~ Mc ) for > (37) 

2cosh(n/j c ) r 
Af c sinh(n Mc ) iU1 M — Mc- 



The large fi> c N c , fiN c limit of the microscopic phase quenched partition function (j34|) for /i < /i c is given by 

16n 4 (/i 2 — /Lt 2 )/x 
whereas the two-flavor partition function in this limit reads 



Z pq ( Mc)M )- 4 r„2_„2 W > (38) 



Z ^ 2(,1 ^ ) = W (39) 

so that the average phase factor simplifies to 

(e 2ie ) pq = l-4 for ^ </ic - ( 4 °) 

This result is in agreement with the mean field result (flU)) . Indeed, the denominator of (|38[) is the product of the 
Goldstone masses fi c — fj,, n c + /i, /i c and /i c , whereas the denominator of (|40[) is the product of four Goldstone masses 
fi c . The free energy of both partition functions is equal to 2nN c fi c . 

For /i > fi c , the large N c limit of the two flavor partition function is still given by (|39|). but the large N c /i limit of 
the microscopic phase quenched partition function is now given by 



& 2nN a p, 

I6n 4 (fi 2 — /i 2 )/i 2 ' 

resulting in the average phase factor 



^pq = TeZJTTTa ..2\..2 ' (^1) 



Pq = %^ = 4(4 - lJe-^-l"-^ for M > Mc • (42) 



2i0\ _ Z ^/=2 _ M /M n -2nJVo( j u-Mo) 



This result can also be obtained from the small n/i, n^t c expansion of the large N c limit of the average phase factor 
(see Eq. (|37|)) at fixed nN c fi c . We conclude that the sign problem becomes exponentially hard for fi > fi c . 

In Fig. [T] we show the average phase factor for different values of n and N c and a critical chemical potential of 
He = 0.1. In the thermodynamic limit at fixed N c , the average phase factor is one for // < /i c and jumps to zero for 
fi > /i c . In the right figure we observe a rapid convergence to the mean field result discussed in section [II CI 



B. Average phase factor for full QCD 

The average phase factor with the ./Vy-flavor quark determinant as weight is defined by 

„2i8\__ _ 1 f jTT detD ^.N 



)N f 



dU- — — -det" f D 



Zjv> (Me A*) J det.Dt 

4-4 r / dU ,) Y % Jdet^(l - Ue^- n ^)det N f(l - u^e~ n ^ n ^). 
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FIG. 1: The average phase factor for one-dimensional QCD with gauge group U(N C ) calculated as the ratio of the two flavor 
flavor partition function and the phase quenched partition function. In the left figure the number of lattice points varies as 
indicated for fixed N c = 3, and in the right figure the number of colors varies as indicated for fixed n = 4. 



We evaluate this integral exactly for Nf = and Nf = 1 but only give its large N c limit for other values of Nf. 
Using the CFZ-formula (fl~5|) [la ] we obtain in the quenched case 

< e )jv /= o- ( 1 _ e -2n ; , e )2 + e (l-e 2n ^)(l-e- 2n ^) M ( } 

The large iV c limit of this result coincides with the large N c limit of the phase quenched average phase factor. In the 
microscopic limit where fi c nN c and finN c remain fixed for N c — > oo we obtain 

2 2 

(e ;jv/=o - 1 2 2 e • ^> 

In the large nN c /j, c limit the second term does not contribute resulting in 

<e 2i % /=0 = l-4, (45) 
K 

in agreement with the results obtained in [4[ (see subsection III C|) . In the thermodynamic limit at fixed -/V c the average 
phase factor converges to one for fj, < fj, c . 
After rewriting the determinants in (|43|) as 

det(l-^t e -n^.) ^ ^ det (l - Ue*M+*Mc) 

det(l - [7t e »M-nMc) ^ det(l - J7e- n ^+"^ c ) ' 1 J 

the denominator of (|43[) can be expanded in powers of £/ for \i > fi c . For Afy = only the constant term in the 
integrand yields a non- vanishing result. We thus find 

(e 2i % /=0 = e- 2nN ^ for ji > fi c , (47) 

so that the average phase factor vanishes in the large nN c /ji limit. 

In Fig. [2] we show the quenched average phase factor for various values of n and N c . Also in this case the average 
phase factor jumps from 1 to at ji = /i c in the thermodynamic limit. The convergence to the mean field result for 
increasing N c (right figure) is very rapid. This can be understood from the expansion of the microscopic result given 
byl-( M 2 /M?)(l-f>M) 2 + --, 
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FIG. 2: The quenched average phase factor for one-dimensional QCD with gauge group U(N C ). In the left figure N c = 3 and 
the number of lattice points is as indicated. In the right figure the number of colors varies whereas n = 4. 



Finally we calculate the average phase factor for an arbitrary number of flavors in the large N c limit. The CFZ 
formula can again be applied after inserting the factor 

det*'(l - Ue- a )det N f(l - ^V' 3 ) 

in the denominator and taking the limit a — > oo and (3 — > oo at the end of the calculation. 

For fi < fi c , the leading order large N c result is given by the identity permutation in the sum over permutations in 
the CFZ formula. This results in 



(T _ g-2n/j c ^2 



Nf+1 



for fx < fx c . 



(48) 



The microscopic limit of this result is given by 



\ e )N f - 



i-!L 



Nf+1 



for fj, < ji c 



(49) 



in agreement with the discussion given in section III CI 

For [i > ji c the CFZ formula can again be applied after rewriting the partition function according to (|46|) and the 
insertion (|48[) . Because of degenerate prefactors, all critical chemical potentials that occur in the combination fi + \i c 
have to be taken different. After carefully taking limits we find 



^jv>+i|i» (Mo Mj 



N f \' 



-2fmN c 



(1 - e - 2 ^ n ) N f 



(1 _ p -2n^N f 



(50) 



The average phase factor given by 



{ e )N f = 



Zn, 



- 2nN ^ for LL> 



(51) 



vanishes in the thermodynamic limit. The exact result becomes rather cumbersome even for small values of Nf. As 
an illustration we give the result for Z 2 \i* (/x Cj m) m Appendix A . 
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V. AVERAGE PHASE FACTOR FOR SU(N C ) 
To calculate integrals over SU(N C ) we use the identity 

oo . 

dU---= V / dUdct p (U) ■ ■ ■ . (52) 
For the partition functions discussed below, with a few exceptions the sum truncates to a small number of terms. 



A. SU(N C ) QCD Partition Function 

Let us first calculate the Id SU(N C ) QCD partition function for Nf = 1. The partition function is defined by 

Z S x {Nc) (n c ,n)= ( dUdet N W = V f df/det p C/e nAr ^'=det^(l + C/e n ^~™^)det^(l + C/ t e^""^). (53) 



>SU(N C ) p^-^JUiNc) 

For p > Nf the integrand can be rewritten in terms of determinant of matrices U only 



Zjv /)P (Mc,m) = / dU det p U det N f D 



U(N C 

dUe nN " N f^det p ~ N f (U)det Nf (1 + Ue n ^- n ^)det Nf (U + e -w-"*>) . (54) 

U(N C ) 

For p > Nf it follows immediately that the integral vanishes. For p = Nf only the constant term inside the 
determinants contributes to integral resulting in 

ZNs,p>Nf(pc,v) = 0. (55) 
For p < we combine det p J7 = det~ p [/t with the factor det^ (1 - Ue nfJ >- n ^). We then find 

Zjsr fiP --N f (fJ'c, fj) — e nJV/JVcM , 

ZN f ,p<-Nf(pc,V) = 0. (56) 

Using that the hermitean conjugate of a unitary matrix is also unitary we obtain 

Z Nf -p(n c ,fi) = Z NflP (iJ, c ,-fx). (57) 

For \p\ < Nf the integral can be calculated by means of the CFZ formula by choosing all [i c different and introducing 
the limit 

det U = lim det(U - e~ a ). (58) 

a — >oo 

The limit of degenerate /i c and a — > oo is taken at the end of the calculation. We will only give exact results for 
Nf = 1 and Nf = 2. 

Using (|55l56p and the result for U (N c ) we obtain for Nf — 1 

ztf^M = e-W + e*W + Bmh(n(jV e + IK) 
f smh(n/i c ) 

in agreement with earlier work by Bilic and Demeterfi [24j]. 

For Nf = 2, the only integrals that have not yet been calculated are those for p = ±1 which are related by (|57[) . 
For the p = 1 contribution we find 



Ziv /= 2, P =i(Mc,/i)=e n2V <^-^ 



(60) 
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For fi < fi c the large nN c fi limit of the Nf = 2 partition function is dominated by the p = term, whereas for 
/i > /i c the large nN c fi c limit is given by the p = — 2 term. For /i = fi c the large iV c limit is dominated by the p = — 1 
with 1/JV C corrections from the p = and p = —2 terms. We thus find as leading large N c result 



, , , g2n(We+2) i Ue 

Z AT f =2 (MoM) = — . , 4/ \ fOT M<Mc 

lb sinn \nfi c ) 



4^2 c) (McM) = e 2 " w ^ for /i> Mc , 

4^2 c) (McM) = 4 C sinh 2 (n/i } [l + O^c" 1 )] for // = M c- (61) 

Both for fi < fi c and /i > /i c the terms that are canceled by the integration over the unitary group are subleading in 
the thermodynamic limit. In other words there is no serious sign problem. 

By inspection one can easily show that the dominance of the p = terms for fi < fi c and the p — —Nf term for 
fi > fi c is a feature of the large iV c limit that is valid for any number of flavors. Therefore, for // < fj, c , the chiral 
condensate is the same as for the U (N c ) theory. For fi > /i c though, the dominant p = —Nf term is mass- independent. 
This results in a vanishing chiral condensate so that the large N c limit of the U(N C ) theory and the SU (N c ) theory 
is different. 



B. The Phase Quenched Partition Function for SU(N C ) 

A second ingredient for the average phase factor is the phase quenched SU(N C ) partition function. Using the same 
arguments as for the full QCD partition function, one easily obtains 

^l+lMp|>2(A t c Mj = 0; 
^1+1*,| p |=2(A*c,M) = 1; 

Zi+i*,-p(a*c,m) = Zi + i* iP (h c ,-ij,). (62) 

The partition function for p = is the U(N C ) partition function which was already given in (|33[) . What remains to 
be calculated is the partition function for p = +1. Using the CFZ formula one easily arrives at 

z _ , j = 1 / sinh(n(iV c + 2)( Af + Mc )) _ sinh(n(iV c + 2)(/i - Mc )) \ 

4sinh(n/^) sinh(n/i c ) \ sinh(n(/i + /i c )) sinh(n(/x — /i c ))) / 

Both for fi < [i c and /i > /i c , the leading large iV c result of the phase quenched partition function resides in the 
p = term given in Eq. (|35|) . 



C. Average Phase Factor for SU(N C 



The average phase factor can be either obtained from the ratio of the full QCD partition function and the phase 
quenched partition function or can be calculated with the full QCD partition function as weight (see Eqs. @ and 
([3])). In the first case we find for large N c , 



sinh(n(/^ c — fi)) sinh(n(/^ c +/^)) 
sinh^ (nfi c ) 

16sinh(n(/i — /i c )) sinh(n(/i c + fi)) sinh 2 (n/i) 



for < He, 

for ^ > jU c , 
for /i = Mc- 



(64) 



The microscopic limit of the average phase factor is obtained by expanding this result for small fj, c , /i. We again find 
the usual mean field result 



a 2i8\SU(N c ) 
- /pq 



6(ji c - fi). 



(65) 



In the thermodynamic limit, the average phase factor converges to one in both for fi < fi c and /i > fi c . 
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FIG. 3: The phase quenched average phase factor for one-dimensional QCD with gauge group SU(N C ). In the left figure 
N c — 3 and the number of lattice points is as indicated. In the right figure the number of colors varies whereas n = 4. 



Exact results for the average phase factor are displayed in Fig. [3] At finite N c , we observe a rapid convergence 
to the asymptotic value of 1. At fixed n, the approach to the large N c limit is slow. Also clearly visible is that for 
fi ~ [i c the corrections to the microscopic limit are large. As was already discussed before, this is due to the p = — 1 
contribution. 

In the quenched case the average phase factor for SU (N c ) defined according to Eq. ([3]) can be written as 

{£ )N f-° ~ p ^Jueu(N c ) 1 j det(l + [/ e -«™c)det(l + t/te^-^ c )- ^ 

The integrals can again be calculated by means of the CFZ formula. For N c > 3 one can easily show that the limit 
introduced in (|58p gives vanishing results for \p\ > 2. However, for N c = 1 and N c = 2 additional terms may contribute 
to this limit. For N c = 2, where the average phase factor is equal to one, contributions for \p\ > 2 vanish for /i < /i c , 
but terms with p < —2 are nonzero for fj, > [i c . The formulae for \p\ < 2 given below are valid for N c — 2. For N c = 1 
the integrals in (|66p only vanish for p > 2 and fi > fi c . Also the general answer for p = — 2 is not correct in this case. 
For p = — 1 we obtain 




i 1 Jc e ' 2 " M) [(1 " e~ 2 "^~ 2 "^) - e - 2nN ^ ( e - 2n ^ - e' 2 "^)] for /i > Mc- 



(67) 



(e 2l % /=0 = ^ (68) 

for [i > fi c . 

Finally, for p — —2 the result is 

{0 for fi < /i c , 

(69) 
e -2n M )(l _ e -M^+Mc))( 1 _ e 2n(Mo-M))(l_ e -2«M) f or M > Mc . 

For p = 2 the average phase factor vanishes. Notice that the relation 

(e 2l6 ) Nf =o,- P (Vc M) = (e 2te )N f =oAVc -M) (70) 
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FIG. 4: Average phase factor for one-dimensional QCD with gauge group SU(N C ). In the left figure N c = 3 and the number 
of lattice points is as indicated. In the right figure the number of colors varies whereas n = 4. 



is also valid in this case. It can be used to derive results for negative values of pL. 

In the large N c limit, the average phase factor is dominated by the p = contribution for fi < /i c and by the p — — 2 
contribution for /i > [i c . For [i = /_t c , the p = 0, p = —1 and p = —2 contributions are of equal order in N c , but only 
the p = — 1 term is nonvanishing. As large N c limit of the average phase factor we thus find 



, p 2ie\SU(N c ) _ 

\ e /N f =0 ~ 



sinh(n(/j c — /^)) sinh(n(/i c +^i)) 
sinh 2 (rtfi c ) 



for (i < fi c , 

16 sinh(n(/i — fi c )) sinh(n(/i c + /i)) sinh 2 (n/x) for fi > fi c , (71) 

1 — e~ 4nf * c for fj, = fi c , 



which is the same expression as obtained in (|64(1 for the phase quenched average phase factor. The microscopic limit 
of the average phase factor is again given by the small /i-/i c expansion of this result which is equal to the usual mean 
field result 

(e->™=(l-|)^ c -M). (72) 

In Fig. 4 we show the exact result for the SU(N C ) average phase factor. We observe a rapid approach to the 
thermodynamic limit at fixed N c just like in Fig. 3. 

In the large N c limit the result for the quenched average phase calculated according to (f6"6"]) is also given by . As 
already was argued in the introduction, we conclude that QCD in one dimension with SU(N C ) as gauge group does 
not have a serious sign problem, not only for /x < fi c , but also for p > fi c . Indeed one-dimensional QCD at nonzero 
chemical potential could be simulated reliably by the Glasgow method [26| ■ 



VI. THE SIGN PROBLEM AT NONZERO TEMPERATURE 



The results of the previous two sections show that the sign problem is severe for QCD in one dimension with 
gauge group U(N C ) and fi > fj, c . Although this theory does not have a critical temperature, it is clear from Fig. 1 
that the average phase factor for /i > fi c increases significantly for higher temperatures (i.e. for lower values of re). 
However, this is a geometric effect due to the Boltzmann factor. Recent lattice simulations [13, El, E3, HB, HI, HH, 
|40| . l4ll . l42l [43l [44], |45|, H^, |4?J suggest that the sign problem becomes much milder for T around T c . Because of the 
absence of a critical temperature, this temperature dependence of the average phase factor cannot be investigated in 
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FIG. 5: Scatter plot of the eigenvalues of the random matrix Dirac operator for =0.3 and T = 0, 0.7, 0.9, 1.0, 1.1 . The 
quark mass at m — 0.1 is indicated by the black dot. 



one dimensional QCD. Instead we study this problem in a schematic random matrix model at nonzero temperature 
and chemical potential [3(| 0, . This model has a severe sign problem as well as a critical temperature and its 
eigenvalues are scattered in a finite domain of the complex plane, just like in QCD. The quenched average phase factor 
in this model is given by the matrix integral 



with Dirac operator given by 



, iC + v + it , 
D - I iCt + n + it o 1 (74j 



and t is the traceless diagonal matrix t = diag(— T, • • • , — T, T, • • • , T). The integral is over the real and imaginary 
parts of the matrix elements of the complex N x N matrices C . The phase quenched average phase factor is defined 
by 



fdCdet 2 (D + m)e~ NTrCC " 



y /pq fdCdet(D + m)det(Di + m)e- Nr * cc1 K ' 



If at all possible, it certainly takes a significant effort to evaluate these integrals analytically. Therefore we have 
studied the quenched (I73|) and phase quenched ((75|) average phase factor numerically. In Fig. 5, we give scatter plots 
of the eigenvalues of the random matrix Dirac operator. For each of the temperatures T = 0, 0.7, 0.9, 1, 1.1, the 
figures show results of 40 matrices with TV = 100 and // = 0.3. The black dot represents the quark mass for which the 
average phase factor is calculated. In Fig. 6 we show the ^-dependence of the phase quenched (left) and quenched 
(right) average phase factor obtained by averaging over 10,000 to 80,000 matrices C with N = 100. The results for 
the phase quenched average phase factor for T = and T = 0.7 and /i > fi c are not displayed because they do not 
converge because of the severity of the sign problem. The average phase factor depends only weakly on the size of 



1G 




the matrices. The conclusion of Figs. 5 and 6 is that the sign problem is not serious when the quark mass is outside 
the domain of the eigenvalues. This confirms the conjecture that was made in 0, [H, H3j- ^ also agrees with the 
observation in [30| that the free energy of the full theory and the phase quenched theory are the same for T = if 
the quark mass is outside the domain of the eigenvalues (this is the case for /i < fi c ). The /^-dependence is due to the 
curvatures at the saddle point. 



VII. COMPARISON WITH THE QCD 3 PARTITION FUNCTION 

According to the general arguments given in section II B, the microscopic limit of the one dimensional QCD 
partition function at nonzero chemical potential is equal to the microscopic limit of the QCD3 partition function at 
zero chemical potential, but with shifted quark masses. In this section we show this explicitly for some of the results 
derived before. We first discuss the microscopic limit of the QCD3 partition function which can be derived from either 
a chiral Lagrangian or a random matrix model with the same global symmetries [28| . For a discussion of QCD3 in 
terms of chiral Lagrangians we refer to [5l[ . 



A. Random Matrix Model 



From the continuum Dirac operator Q it follows that the random matrix model that describes the fluctuations of 
the low-lying Dirac eigenvalues is given by 

Z = J dHP(H) det(iH + [i + m) det(iH -m + (76) 

where the probability distribution can conveniently be taken to be the Gaussian distribution 

P(H) = e ~^TrH\ (7?) 

Such random matrix partition functions have been studied elaborately in the literature to analyze the microscopic 
limit of QCD in three dimensions [H, [52|, HH, [U, HH . The microscopic limit of the partition function (|T6[) with quark 
masses given by M — diag(— mi, • • • , —rriN f , mi, ■ ■ ■ , mjv,) is equal to [53|, HH 

ZQCD3 ~A(M) det {A(-m) A(m) J' (78) 
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where the matrix elements of the Nf x Nf matrix A(m) are given by 

A{m)ki 

and A is the Vandermonde determinant 

2N f 



m l k e- mk 



A(M) = H(M kk -M u ). 



(79) 



(80) 



k>l 



Although general expressions for QCD3 partition functions with an arbitrary number of bosonic and fermionic 
determinants are also known [56| , for our purposes we only need the partition function 



J N f +2\2 



n f , ;/r^ s , ;/r^ .det(D + z 1 )det(D-z 2 
det(£ + m k ) det(D - 



\fe=i 



--2 



(81) 



This partition function was evaluated in [57j by means of the supersymmetric method and is given by 



JV, 



Z 



N f +2\2 



(z\ - zi){z 2 - z 2 ) exp7V c (zi + z 2 ) -q (m k - z x )(m k - z 2 ) ^ 



N c (z! + z 2 ) 1J- {mk - Zi)(mk - z 2 ) 



sinh N c (rrik +rrii ) sinh N c {m k +zi ) 

N c {m k +mi) N c (m k +zi) 

N c {z2 J rrai) sinh A^ c (21+^2) 

iV c (^2+m,) iV c («i+z 2 ) 



-Af c (zi-Zl) + AT c (z 2 -22) 



■det 



sinh iV c (771 £ + mi) 



N c (m k + mi) 



(82) 



The last term is a so-called Efetov-Wegner term, and the Vandermonde determinant is over positive masses only 

N f 



A (m) = JJCm/s - to/). 



(83) 



k>l 



B. Microscopic Limit 



We will now show that the microscopic limit of the partition function (|25|) is equal to the microscopic limit of the 
QCD3 partition function. To this end we multiply row k of the first Nf rows of the determinant by exp(— N c m- fc/2) 
and row k of the second Nf rows of the determinant in (|25|) by exp(— N c m +k /2). In the microscopic limit we keep 
N c m±^ k fixed so that we can expand the masses that do not occur in this combination. By subtracting successive 
columns starting with the first one we obtain 



Z Nf = 



1 



III 



<k<l<2N f 
-jVem-i/2 



(M, - M k ) 



(84) 



m_ ie 



-7V c m_ i/2 



m N f -l e - Nem _ l/2 



, N f +N c m^ 1 /2 



-N c m- N ,/2 -N c m- N ,/2 
e /' m^N f e f 



-N c m + i/2 



-N, 



m + JV/ /2 



777+ ie 



-JV C 



-1/2 



m »f-yN.rn- Nf /2 m ^ e +iV c m_ N/ /2 

777+ / 1 Ll e- JV = m + 1 / 2 m^e 



v / +iV c m + i/2 
■ l r 



-A', 



c m + Nf j2 



Nf-1 -N a m +N ,/2 N f +N c m +Nf /2 

m +Nf e + "f rn + ' Nf e +"f 



2AT/-1 w c m-i/2 



2JV/-1 N c m^ N j2 
m 2N^-l eNaTn+l/2 



2N S -1 N c m +Nr /2 
m +N f e * 



where the masses Mk arc defined below Eq. ([23)) . By multiplying the columns / = Nf + 1, • • • , 2Nf by (— l)'^ 1 and 
introducing microscopic masses, we obtain exactly the expression for the QCD3 partition function given in (|T8[) . 

For masses — (/j c + /i), — (fx c — )i), fJ, c — fi, fi c + fi corresponding to the microscopic limit of the two-flavor phase 
quenched partition function given in (|34[) the determinant in (|84[) is given by 



477 M 2 (e 2 " JV ^ + e -2nN^ c) _ 4n/i 2 (e 2„2V oM + e -2nN c ^ + g^a _ ^ 

and the prefactor is equal to 

J] (Mi - M k ) = 64n 6 (/^ - at 2 )^ 2 . 

l<k<l<2N s 



(85) 
(86) 
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FIG. 7: The number variance for n — 4 and N c = 100. 



Their ratio coincides with the microscopic phase quenched partition function given in (|34[) . 

As second example, we consider the quenched average phase factor for U {N c ) given in (|4"5|) . For Nt — the partition 
function (|82|) simplifies to 

. e N ^ Zl+ ^ sinhNJzi + z 2 ) M t, * , + * \ , , 

Z 2 \2 = (zx - zi)(z2 - z 2 ) c{ ' V + e -^(*i-*i-«+«). 87 

N c (zi + z 2 ) N c (zi + z 2 ) 

The masses in this partition function corresponding to (|43p are given by 

z\ = nfx c + rifjb, z 2 — nfj, c — n/i, 

z 2 = n^L — njjL c , z\ — —n[L — n\i c . (88) 

Substituting them into ([57)) we obtain 

Z 2 \ 2 = 1 - 4 + 4 e " 4JVcnA1 (89) 
which is exactly the microscopic limit of the quenched average phase factor (|4"3"|) given in Eq. (|4"4f . 



VIII. DIRAC SPECTRUM 

We consider the Dirac operator (|4]) in a gauge where all gauge fields except U n \ = U and Ux tTl = W are equal to 
unity. If the eigenvalues of U are equal to expiiOk) the eigenvalues of the Dirac operator are given by 

I 2iri(fc + l/2) + i£), 2wi(k+l/2)—iB, 

X k ,i = -(e s +^-e s f), fc = i,..., n , l = l,---,iV c . (90) 

Contrary to QCD at /i =/= in more than one dimension, where the eigenvalues are scattered in the complex plane 
[27l I30I [58l . [59l . [60l [6l| , the eigenvalues are located on an ellipse in the complex plane with real and imaginary parts 
related by 



2 / T _/\ \ \ 2 

I 



/ Re(A fc ,;) y / Im(A fc , 



= 1. (91) 
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A. 



Universal Fluctuations 



In section III Bl we have shown that the continuum limit of the staggered Dirac operator in one dimension is in the 
same universality class as QCD in three dimensions. The continuum limit of the staggered three dimensional Dirac 
operator, though, is in the chiral symmetry universality class of QCD in four dimensions. For example, in [62| this 
was found for the distribution of the small Dirac eigenvalues. The reason for the chiral structure is that the staggered 
Dirac operator only couples even and odd lattice sites. This is also the case in one dimension, but the off-diagonal 
blocks, containing the gauge fields, are the same in the continuum limit (or occur in the combination U n + U^ n ), 
resulting in a two flavor theory with opposite masses (see ([12")) ). The secular equation is given by 



This corresponds to the superposition of the spectrum of do + Aq + \i and — do — iAo — fx. In the domain of the Dirac 
spectrum where do can be neglected, the eigenvalues on each of the lines ±fi are given by the Hermitean random 
matrix ensemble Aq. For fi = we have the superposition of the Hermitean random matrix ensembles of Aq and — Aq. 
From the eigenvalues (|90|) it is clear that the staggered lattice Dirac eigenvalues show a similar superposition. One 
ensemble is localized on the left half of the ellipse and the other ensemble on its right half. Superpostions of more 
ensembles arise for correlations on scales larger than N c level spacings. 

Eigenvalues of U (N c ) matrices are correlated according to the circular unitary ensemble. Therefore, the eigenvalues 
Ak,i are correlated according the Gaussian Unitary Ensemble (GUE) on a scale where the variations in the average 
spectral density can be neglected, i.e. on a scale much less than N c level spacings. By correcting for the variation of 
the average level spacing, the scale on which universal random matrix correlations are found can be extended, but this 
scale should always remain well below N c level spacings. Then we are necessarily within one of the n successive copies 
of the eigenvalues along the ellipse. At a scale of a finite number of average level spacings, we expect convergence to 
universal random matrix correlations in the large N c limit. 

As measure of the correlations between eigenvalues, we use the number variance S 2 (n) defined as the variance on 
the number of eigenvalues in an interval that contains n eigenvalues on average. We calculate the number variance for 
intervals starting from zero. Because of the superposition of an enemble and its negative, this is equal to the number 
variance of a single ensemble for an interval that is symmetric about zero which is given by the GUE result. In Fig. 
Owe show the number variance obtained from an average over 10 5 gauge field configurations for N c = 100 and ^ = 0. 
The number variance is calculated for an interval starting at zero. The discrepancy between the analytical result and 
the GUE (solid curve) is only barely visible. 

The increase of the domain of validity of chiral random matrix theory with increasing N c was also observed for 
Dirac operator of QCD in four dimensions [63| . This can be explained as follows: The number of eigenvalues that is 
described by random matrix theory is equal to ~ F%y/V (with V the space-time volume), but F% scales with N c in 
the large N c limit. 



We first discuss the case of fi = 0. Then the average spacing of the eigenvalues scales as l/(nN c ). This means 
that the chiral condensate develops a discontinuity at to = when nN c — ► oo. Two cases of interest where a nonzero 
chiral condensate can be obtained are the limit n — > oo for fixed N c , i.e. at zero temperature, or at fixed n in the 
limit N c — > oo. In both cases the chiral condensate is discontinuous for nN c — > oo when the quark mass crosses the 
imaginary axis where the eigenvalues are located. However, Goldstone bosons and universal behavior is only found 
in the large N c limit. At fixed N c , chiral symmetry breaking is not associated with universal behavior. 

For /i ^ the Dirac eigenvalues are located on the ellipse (|9"Tj) . When the quark mass is inside the ellipse of 
eigenvalues the chiral condensate is zero for nN c — > oo for the SU(N C ) partition function and the phase quenched 
U(N C ) partition function (24)]. For gauge group U(N C ) with Nj > 1 flavors, though, the partition function is //- 
independent so that the chiral condensate is also nonzero if the mass is inside the eigenvalue ellipse. The Banks-Casher 
formula fails in this case which is known as the "Silver Blaze Problem" . The resolution is the same as for QCD 
in four dimensions [3lj : the re- weighted eigenvalue distribution defined as 



det[(<9 + iAo + M + A )(^o + iA a + fi - A)] = 0. 



(92) 



B. Chiral Symmetry Breaking at fi ^ 




fc 



(93) 



shows oscillations with an amplitude that diverges exponentially with n and a period that is proportional to 1/ra. 
Below we will illustrate this for U(l). 
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For N c = 1 one easily derives that the spectral density is given by 

a nfj, c _j_ p — rifle _ f>n(icc+fi) _ g — n{ia+n) 



p tull (z) = 4 7 : ^ = T . —8( r -n), (94) 

where z is parameterized as 



- -{ e r+la - e- r ~ la ), r > 0, a € [0, 2tt], (95) 



2 



and we have used that Z^ f —\ = 2cosh(n/^ c ) (see eq. ((5)). For p > /i c this spectral density has oscillations with 
an amplitude that diverges exponentially with n and a period of order 1/n which are the essential properties of the 
spectral density of the Dirac operator of QCD in four dimensions |3ll [6a. 66]. It can be decomposed as 

p iM (z)=p«(z)+p° sc (z). (96) 

The average quenched spectral density is defined by 

P%z) = Y [ d9-j28 2 {z~\ k ) (97) 
71 J n fe=i 

with the eigenvalues Xk given by Eq. ([90]) . After changing variables according to (|95[) and integrating over 8 we find 

P<i{z) = -k {e*- + e-ir + e ^ + e-^) 5{r ^ (98) 

The oscillatory part of the spectral density, p osc (z) is equal to the difference p full (z) — p q (z). 
The chiral condensate is given by 

and can also be decomposed as 

E full (m) = S q (m) + E osc (m). (100) 
The Jacobian for the transformation of d 2 z to drda is given by 

d 2 z = i(e 2r + e~ 2r + e 2la + e - 2la )drda, (101) 
so that the chiral condensate after integration over r can be simplified to 

.full/ 1 ( 2V , e n/ic + e ~ nAlc - e"( M+ ^ - e - n t ia+ ^ _tanh(np c ) 



S lull (m) - - / da-r — rr—. - ^ = ir^p^. (102) 

v ' 7T J ( e ^e -i- e -^c)(eC+« - e -c- M - 2m) cosh/x c 

This result remains finite for n — > oo. We remind the reader that the mass is parameterized as m = sinh/i c . This 
expression represents the resolvent at the quark mass which was evaluated in [49j and is in agreement with earlier 
work d, [24[ . Decomposing the spectral density and the chiral condensate according to we obtain 

E osc (m) = <9(sinh(/i) - m) £ full (m), 
E q (m) = 0(m- sinh(^))E full (m), (103) 

so that when the quark mass is inside the ellipse of eigenvalues, the entire chiral condensate is due to the oscillatory 
part of the spectral density. The discontinuity in the chiral condensate is reminiscent to a Stokes phenomenon. The 
alternative to the Banks-Casher relation proposed in [3l[ is also at work for QCD in one dimension. This solves the 
"Silver Blaze Problem" [641 ]. 
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IX. CONCLUSIONS 



We have studied QCD in one dimension at nonzero chemical potential. Both the full theory, its quenched and 
phase quenched versions and gauge groups U(N C ) and SU(N C ) have been considered. In one dimension, the QCD or 
QCD-like partition functions can be reduced to a single matrix integral, which because of recent advances by Conrey, 
Farmer and Zirnbauer, could be evaluated analytically. In this paper we have analyzed the small mass behavior of the 
partition function, the nature of the sign problem, and the relation between the Dirac spectrum and chiral symmetry 
breaking. 

To put our results in perspective, we emphasize that QCD in one dimension, is quite different from QCD in 
more dimensions. In particular, we wish to mention the following three points. First, instead of being scattered in 
the complex plane, the eigenvalues of the staggered Dirac operator are located on an ellipse in the complex plane. 
Second, phase transitions can only take place for zero temperature or for N c —>oc. At the critical chemical potential, 
a transition from the vacuum state to a state of free quarks takes place. Third, color singlets are made out of 
noninteracting quarks so that the critical chemical potential is given by the quark mass. In particular, the critical 
chemical potential for the meson state and the baryon state is the same. 

Because the Dirac eigenvalues are located on a curve, the large N c limit of staggered lattice QCD in one dimension 
is in the same chiral symmetry class as QCD in three dimensions. We have shown this both by an explicit evaluation 
of the partition function, and by analyzing the fluctuations of the Dirac eigenvalues for large N c . We have used this 
equivalence to explain the behavior of the partition functions and the average phase factor in the microscopic limit. 

Contrary to QCD in four dimensions, the sign problem for QCD in one dimension is not severe in the case of gauge 
group SU(N C ). One reason is that the critical chemical potential for full QCD and phase quenched QCD is the same, 
so that the parameter domain m n /2 < fj, < mjv/3, where the sign problem becomes severe in four dimensions, is 
absent in one dimension. A second reason is that for /i > fi c both the full theory and the phase quenched theory 
become a theory of free quarks with the same free energy in the thermodynamic limit. 

For gauge group U(N C ), on the other hand, the sign problem is severe when /i > both in one dimension and in 

four dimensions. The reason is that the U (N c ) theory does not have charged excitations, whereas the phase quenched 
theory has charged mesons resulting in a phase transition at (J, = rn^/2. We have evaluated the average phase factor 
both by averaging with respect to the phase quenched partition function and the full partition function, and similar 
conclusions have been reached. 

The condititon \i > for having a severe sign problem can be rephrased as the quark mass being inside the 

ellipse of eigenvalues. In more dimensions this condition is that the quark mass is inside the support of the Dirac 
spectrum. It also applies to nonzero temperature as we have demonstrated explicitly in the framework of a chiral 
random matrix model. 

Also for U(N C ) QCD in one dimension, the chiral condensate is discontinuous across the imaginary axis in spite 
of the fact that there are no Dirac eigenvalues. This can only mean that the phase of the fermion determinant 
is responsible for the discontinuity. This is what happens in four dimensions where the discontinuity is due to a 
contribution to the spectral density that oscillates with a period of the inverse volume and amplitude that diverges 
exponentially with the volume. Exactly the same mechanism, where the discontinuity in the chiral condensate arises 
due to a Stokes like phenomenon, is at work in one dimension. This suggests that this is a universal mechanism for 
theories with a sign problem. 

We end by repeating that QCD in one dimension with SU (N c ) as gauge group has no serious sign problem. Our 
hope is that part of this conclusion translates to four dimensions ameliorating the sign problem for fi > M/v/3. 
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APPENDIX A: AVERAGE U(N C ) PHASE FACTOR FOR N f = 1 

In this appendix we give explicit results for the average phase factor for Nf = 1 defined by 



2,i) Z 2 \\* (fi c , aO 



The partition function ^jv*=i(Mc) A 4 ) was already given in Eq. 1(8"]). The numerator can be obtained from the CFZ 
formula. After taking the limit of degenerate critical chemical potentials at the end of the calculation, we obtain for 
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^2|i.(McM) = e-'^^e- 4 "" 1 , J 1 6 „ L + e -3n(JV c +iK if ? > 

+e -„(iV +l)^ c [/q + ^ + ^2^ (A2) 

with 

2g2n(/i+/i c ) _|_ 2g2n(/i-/id) _ g4riM _ g-4n/j , ^g2n^i _|_ ^g-2n/j _ g4n(/i— /j c ) _ e 4ra(/i+/j c ) 2g6njll c — 2n/4 _|_ 2g6n/l c +2n/l 



/o 



/l 



-10e 2 ™^ + 6e 4r ^<= - 2 - 5e 6n ^ - e" 2 ™^ 



(1 _ e _2n ^) 2 (l — e n ^ +Mc ))(l — e n( ^ M ~ Mc )) 

^ (g'T'Mc g— n/x c ^3 " (^^) 

For \x > \x c a similar calculation results in 

2|1 *" Cg (l_g-2n^)2 ^ + JV/ + ° (l_ e 2n^)2 U+ V c j 

(A4) 

with 

1 1 2__ 

^ g2ri/j _ Y p2ri(p+fj, c ) ^ g2n^i c ^ ' 

11 2 

'92 — e 2n(jl-fl e ) — \ e 2nfi _ I g-2n/j c _ ]_ ' ' ' 

We have checked that the results in this appendix agree with the large N c limit of the average phase factor given in 
section IV B. 
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